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Abstract 

The recent introduced dressed coordinates are studied in the path-integral approach. These coordinates are 
defined in the context of a harmonic oscillator linearly coupled to massless scalar field and, it is shown that in 
this model the dressed coordinates appear as a coordinate transformation preserving the path-integral functional 
measure. The analysis also generalizes the sum rules established in a previous work. 

1 Introduction 

In recent works, the concept of dressed coordinates and dressed states have been introduced in the context of 
a harmonic oscillator linearly coupled to a massless scalar field [TJ [21 OH H] . As emphasized in these references, 
the introduction of the dressed (or renormalized [H [5]) coordinate and state concepts is necessary in order to 
give physical consistence to the oscillator-field system as a toy model to describe an atom in interaction with the 
electromagnetic field, where the atom is roughly modeled by the harmonic oscillator. Also, as early stressed the 
introduction of dressed coordinates and states is twofold advantageous. From the physical view point, the dressed 
states behave as the one expected for the physically measurable states: excited atomic states are unstable whereas 
the atom in their ground state and no field quanta is stable. On the other hand, it allows exact computations 
for the probability amplitudes of the different radiation processes of the atom. Indeed, when the calculation is 
performed for weak coupling constant we obtain, for the spontaneous decay of the first excited state of the atom, 
the long know result: e _rt [6]. Furthermore, when applied to a confined atom, approximated by the harmonic 
oscillator, in a spherical cavity of sufficiently small diameter the method accounts for the experimentally observed 
inhibition of the decaying processes [TJ [5] . Besides that, in Refs. [HIE] the extension of the dressed coordinate and 
state concepts for non linear system have been addressed. 

Nevertheless, in all previous works the approach used has been via the operatorial formalism of Quantum 
Mechanics. The aim of this paper is to develop a path-integral approach to the problem. We hope that the path 
integral approach will be more adequate in dealing with the problem of computing the reduced density matrix 
or to obtain the master equation for the atom as is the case in Caldeira-Legget type models [THl HB [H| ■ As in 
previous works, because its exact integrability, the dressed coordinates were introduced in the context of a harmonic 
oscillator coupled linearly to a massless scalar field. 

The paper is organized as follows: In section 2 we introduce the model and compute the propagator by exact 
diagonalization. Section 3 is devoted to establish the dressed coordinates. In section 4 we compute the probabilities 
associated to some physical processes and extend the sum rules found in [3]. Finally, in section 5 we give our 
concluding remarks. Through this paper we use natural units c = K = 1. 

2 The model and its exact diagonalization 

We consider as a toy model of an atom-electromagnetic field system the system composed by a harmonic oscillator 
(the atom) coupled to a massless scalar field. By considering the dipole approximation and expanding in the field 
modes we get the following Hamiltonian [1] 
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H =\{pI + ^lll) + \Y,{PI + "fal) ~QoJ2 C k q k + ±J2^Lq 2 , (1) 

k=l fc=l fc=l k 

where qo is the oscillator coordinate and q k are the field modes with k = 1, 2, ujk — 2ir/L, c k = r/usk, V = ^2gAuj, 
Auj = ujk+i — uJk = 2n/L. With g being a frequency dimensional coupling constant and L the diameter of the 
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sphere in which we confine the oscillator-field system. In Eq. ([I]) the limit N — > oo is to be understood. The last 
term in Eq. (|TJ) can be seen as a frequency renormalization 11 and, it guarantees a positive-defined Hamiltonian. 
Due to the Hamiltonian (fTJ) is quadratic in the momenta and there are not constraints we can write the propagator 
for the system as being 



where q = (qo,qi, 
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and A is a symmetric matrix whose components are given by 
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To compute the propagator, Eq. ©, we introduce the coordinate transformation 
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where is an orthogonal matrix that diagonalize A, 

d = t t at = diag (ng, nl n 2 2 , ■ ■ ■ , nl^nj,) . 

It is easy to show P3] that the eigenvalues of A, £l r are obtained by solving the equation 
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It has shown that such equation has definite positive frequencies ft 2 as solutions PQ. The solutions of the char- 
acteristic equation 1(5)) were used to describe radiation process in small cavities [3J [3] in good agreement with the 
experiment. In [3J this system is also used to describe a Brownian particle coupled to an ohmic environment. 
Replacing Eq. (J6]) in Eq. ([2]) we get 
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Note that in Eq. ^ the functional measure is maintained, this because detT = 1. By using the known result for 
the propagator of a harmonic oscillator we get for Eq. © , 
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1 The components of matrix T and some of its properties are given in the Appendix. As T diagonalizcs the matrix A, the coordinates 
} r are known as normal coordinates 



The spectral function is defined as being 

Y(t) = I dq d qi ■■■dq N K(q, t; q, 0) (f 1) 



which is easily computed expressing the integral in normal coordinates {Q r }, thus we obtain 
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with the energy spectrum being given by 
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In the {Qr } coordinates, the ground state wave function is computed from the propagator given by Eq. (| 10|) by 
taking the limit t — > — ioo 
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thus the wave function of the ground state is 
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We can write the ground state eigenfunction in the original coordinates q^, by using the third equation of ^ 

3 The dressed coordinates 

We have observed the vacuum wave function in the {q^} and normal {Q r } coordinates and we make a question: 
Is it possible to find some new set of coordinates {q^} relate to them what allow us to describe the oscillators 
with their non interacting characteristics? The answer is yes. In such one new set of coordinates the vacuum wave 
function must be given as 

V>oo..o(go, £jv) = cte. exp ^-^ w M {%f^j (17) 

and we named the set of coordinates {q^} as dressed coordinates^ which describe the oscillators as being non 
interacting. In a physical situation we can imagine the atom interacting with a bath of field modes, i.e., an 
electromagnetic field. The experience tell us that the atom does not change his energy spectrum (the energy 
spectrum when it is isolated) and only transitions between the energy levels are observed (absorption and emission 
process). 

Thus we would have for the vacuum state wave functions the following relation 
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First we will look for the matrix transformation between the normal {Q^} and dressed coordinates {f?^}, thus 
we set 
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the quadratic form (| 18[) must be invariant under the linear transformation (|19p , then, to preserve the quadric form 
we set 
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2 The dressed coordinates are a new type of collective coordinates, which are defined from the normal coordinates, that allows a 
correct description for the absorption and radiation phenomena for a given physical system. 



and to achieve the matrix M we use the orthonormal matrix T such that 
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we can show that it satisfies the equation ([20)) by using the orthonormality condition of the matrix T. The 
determinant of the matrix M is shown to be 1, 
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thus the transformation (fT5|) preserve the path-integral measure defining the propagator ^ of the system. And 
the inverse transformation {Q s — > q^} is easily shown to be 
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Then, we come back to the functional integral defining the propagator §§§ in terms of the normal coordinates 
and using the transformation (|22[) we can write 
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from which we can see the Lagrangian in terms of dressed coordinates is given by 
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where we have defined the matrices Z^ v and C M „ 
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such that 

They play the role of renormalization constants and the coordinates {q^} are the renormalized coordinates such as 
it happened in a renormalized field theory. 

To construct the dressed Hamiltonian we first define the dressed momentum p„ canonically conjugate to the 
dressed coordinate q^, thus 
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from which we get 
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Thus, the dressed Hamiltonian Hd is computed to be 
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From (TlTJ|) we write the propagator in dressed coordinates 
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It is simple to show that the spectral function computed in dressed coordinates is the same computed in the 
equation (|12p , thus, the energy spectrum (|13p remains invariant, as it was expected. 
The ground state reads as 
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The Hamiltonian operator in dressed coordinates is expressed as 
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with the coefficients C^. v are given in (f25|) . 
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4 Computing the transition probabilities 

In this section we show what to compute the transition amplitudes of the system using the exact dressed propagator 
(|30p . thus, we are interested in the following quantities, 
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that represents the probability amplitude of the system, initially prepared in the state |mo, mi, mjv), to be found 
at time t in the state \no, ni, njy)- 

Eq. ([33|) can be written in terms of the propagator as 
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Using Eq. (|55|) we can write Eq. (f3"4"| as 
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First we compute A - ^ \- {t). Substituting Vo...On„0...o(x'(x)), #"(x,*;£,0) and Vo...0nv0...o(£'(£)) we have 
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Performing the Gaussian integrals in Eq. (|38p we obtain 
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Using the identity 
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(that holds because {t^} is an orthogonal matrix) where +1 = (lo + h + ... + In) and replacing Eq. (f39|) in Eq. (J3T] 
we get 
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where 
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If instead of integrating over coordinates £ in Eq. (|37[) we first integrate over coordinates x we would get an 
expression similar to the one given in Eq. (|4ip but with Ii rSr replaced with 1[ : 
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Then, since the final result must not depend of the order in which we perform the integrations we must have 
Il rSr = I[ s i and from Eqs. (|42)l and |43|) we conclude that Jj rSr = I Sr i r . 

To perform the integral given in Eq. (|42p we have to use the following theorem 
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Note that the expression in brackets in Eq. (|42[) is a polynomial of degree s r in £ r . Now, if l r > s r , then by using 
theorem (|4"4")l. we get J; rSr = 0. Because 7; rSr = J Sr j r we also get a vanishing result for l r < s r . Then, the only non 
vanishing result is obtained for l r — s r . Using again theorem (|4"4")l we note that the only non vanishing term of the 
polynomial in brackets is the one of highest power. Since the highest power of H n {x) is given by 2 n x n we have for 
Eq. (02]) 
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Using Eq. (gSJ in Eq. |gT} we get 
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where in passing to the last line we have used the identity 
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It is straightforward to establish the following identity: 
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The proof of the above identity follows trivially by using the orthonormality property of the matrix {t r ^}. Writing 
Eq. ([50)1 for indexes and k we have 
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The physical interpretation for the equations above is given as it follows. Let the initial state of the system given 
by \n, lfcj, lfc 2 ...), the atom in the n-th excited level and field quanta of frequencies 0Jk t , ^fc 2 , etc. The probability 
of this initial states to be found in a measurement performed at time t in the state \m, , lfc'-..) is denoted by 

V n . lk \ k 2 (t). We know that 7 ? 1 .' (i) = |/oo(i)| i s the probability of the oscillator to remain in the first excited 

level and P°.Q lfc (i) = |/ofc(i)| 2 is the probability of the oscillator to decay from the first excited level to the ground 
state by emission of a field quanta of frequency Uk- Obviously in this case we have 
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that is nothing but Eq. ([ST) . Also Eq. (1521) can be written as 
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where 
and 

With these identifications the physical meaning of Eq. ([51]) is clear: if initially we have a photon of frequency 
Wfej and the oscillator is in its ground state, then at time t, either the oscillator can go to its first excited level by 
absorption of the initial photon or can remain in its ground state scattering the initial photon to other photon of 
arbitrary frequency. 

Note that in establishing the identities (|53[) and ([5l)l it is used only the orthogonality property of the matrix 
Then, it is natural to ask whether it is possible to compute other probabilities without doing a direct 
computation as performed in last section. The answer is yes. For example, if initially the oscillator is in its second 
excited level and there are no photons, at time t it can happening that the oscillator continues in their second 
excited level, it can go to their first excited level by emission of photon of arbitrary frequency or it can decay 
to their ground state by emission of two photons of arbitrary frequencies ujk 1 and ujk 2 . The respective probabilities 
are denoted by V 2 '-o(t), ^-o" 1 (*) an d ^-o* 1 * a (*)• Obviously we must have 
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Taking the square of Eq. ([55)1 we find 
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Identifying Eqs. ([57]) and ([58]) we obtain 
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As a second example we consider the oscillator is in its first excited state and there is one photon of frequency 
U)ki ■ At time t it can happen that: the oscillator go to its second excited level by absorbing the initial photon; or the 
oscillator remains in its first excited state and the initial photon is scattered to other photon of arbitrary frequency 
ujk 2 ; or maybe the oscillator can be decay to its ground state by emission of a photon of arbitrary frequency u>k 2 
and the initial photon is scattered to other photon of frequency u>k 3 ■ The respective probabilities are denoted by 
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Taking Eq. (J33J) times Eq. jM]) we have 
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From Eqs. (j6"2"|) and (f6"3")) we have 
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And so we can give all the probabilities associated to any decay or absorption process placing in the system. 



5 Conclusions 

In the present work is shown using the path-integral formalism that the dressed coordinates appear as a coordinate 
transformation preserving the quadric form that defines the ground state wave function of the system guarantee- 
ing the vacuum stability and, they also leave invariant the functional measure of the path-integral. Within the 
Hamiltonian formalism it can be shown that such linear transformation, which defines the dressed coordinates, also 
leaves invariant the canonical form of the action. Thus, the dressed coordinates can be also defined via a canonical 
transformation. 

The calculus of the transition amplitudes has been performed using the dressed propagator, being obtained 
the basic formula which defines the sum rules presented in [3], then, the rules have been extended to describe 
other physical processes. In spite of the computation seems very difficult, the dressed coordinates allow to use the 
properties of the Hcrmitc polynomials simplifying greatly the calculus. 

On the other hand, it has been made an extensive use of the model given by Eq. (Q]) to study different physical 
situations, such as the quantum Brownian motion, decoherence and other related problems in quantum optics. In 
such context, it is interesting the computation of the reduced matrix density for the model |T]) in the framework 
of dressed coordinates; the advances in such direction will be reported elsewhere. 

Acknowledgement 

GFH (grant 02/09951-3) and RC (grant 01/12611-7) thank to FAPESP for full support. BMP thanks CNPq and 
FAPESP (grant 02/00222-9) for partial support. 



A The orthonormal matrix T= [t*] 

Because the orthogonal character of the T-matrix its components satisfy 
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and from the Lagrangian ([3]) expressed in terms of the normal {Q r } we get other important relation 
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where u)q have been defined in ((4]) and lo\ = uj\ . Using the equations above we can show the following sum 
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and also compute the elements of the T-matrix 
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